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0\ • Abstract 

\o 

O ! QCD-based analysis of nonfactorizable parts of weak nonleptonic amplitudes is re- 

I~J ■ ported. Nonperturbative effects due to soft gluon exchange play a key role leading to 

Q"^ ■ the emergence of a dynamical rule of discarding 1/Nc corrections. 

^. 1 Introduction 

i-^ ' Factorization is used in nonleptonic decays from early sixties. However, as our knowl- 

edge of QCD and weak decay phenomenology deepens, the simple idea that one must 

^ . factorize two V-A currents composing the effective weak hamiltonian evolves towards 

a rather sophisticated scheme with different ingredients. The purpose of this talk is 
to review recent progress in calculating deviations from naive factorization. We shall 
concentrate here on exclusive decays. ( Inclusive decays are discussed elsewhere.) We 
shall show that the rule of discarding 1/A^c 113' 01 ^^^ a dynamical origin, and is due 
to nonperturbative QCD effects. The key role is played by soft chromomagnetic gluon 
exchange. The resulting picture 0, |, ^, 0, |^ is rather versatile— not all transitions 
are alike in this respect. QCD effects lead to deviations from the naive factorization, 
specific for each channels. These deviations can be estimated in a model-independent 
way. In some channels the situation is close to the predictions of the 1/A^c rule , 
in others — to naive factorization. The degree of cancellation of the naive 1/Nc- 
suppressed amplitudes is different for each channel, so we can call our approach a 
dynamical rule of discarding l/N^. 
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2 The method 

We shall describe our method using the decay B^ — > D^n^ as an example. The reader 
is referred to ref. [Q for details. We are interested in the transitions induced by the 
color-octet times color-octet part of the weak Lagrangian C ~ (cF'^t^fe) {dV^f^u) . In 
order to calculate the transition amplitude, we start with the correlator 



^/5 = I d^x< D\T{C{x),A'^}\B > e^"?^ 



where the axial current A^ = u'-f^'-f^d annihilates the pion. Two key steps are made 
in order to calculate the latter correlator: we continue A^ into the Euclidean region 
—q"^ ~ Q^ ~ 1 GeV^, and write in this region (borrowing some ideas from the QCD 
sum rule method 0) a sum rule for the amplitude M^.i. governed by C: 

A\Q') = M,,J-^ + ... (2) 

where +... denotes the contribution of higher resonances produced by the axial cur- 
rent. Second, we calculate A^ using the Operator Product Expansion. We immedi- 
ately obtain 

A^m = -2z-^^ < D\cTH'^gG:,^b\B > +... (3) 

Svr^ q^ ^ 

where +... denotes higher order power corrections and we retained only the kinemat- 
ical structure proportional to q^ . Comparing the latter two equations and neglecting 
the dots we immediately obtain that the ratio of the nonfactorizable and l/N^. naive 
factorizable parts of the amplitude is 



4vrVr 



(4) 



Note the key distinction from the standard QCD sum rule method: the matrix el- 
ements are taken between hadronic states, not between vacuum states. Using the 
methods of HQET ||T^ it is easy to get for m^^j = \{rn\, — m\). 

It is instructive to emphasize the assumptions and approximations made in eq. 
(4). First, we neglected the corrections due to operators with higher dimensions, and 
contamination with higher resonances. Strictly speaking it is necessary to check that 
the corresponding window exists. This has not been done yet, although arguments in 
favor of smallness of the above corrections in a large class of transitions were given 
in ref. 0]. Second, and this is also important, we started from the theoretical limit 
where Mb — Mo <^ \{Mb + Md). Only in this limit the expansion in eq. (3) goes in 
dimensions, not twists. Moreover, in a number of cases the hadronic matrix element 
in eq. (3) reduces to the known quantity in this limit. Otherwise we would have to 
introduce an unknown function of recoil. Logarithmic corrections due to anomalous 



dimensions are also not included so far (although in the transitions considered in ref. 
[^] they seem to be unimportant). We refer to ref. [Q for the detailed discussion of 
the method and expected uncertainties. The expected accuracy for this particular 
channel is of order one. 

Keeping in mind all these uncertainties— a vast field for future work— one can try 
to extend the method to other weak hadronic decays in a straightforward way. If the 
particle that splits away is not a pion, we do not get a simple l/Q^ term in the OPE, 
but rather a more complicated function. Moreover, the higher power corrections can 
become more important. For example, for the B — ^ DD decays we get the function 
F{Q'^) = l/Q"^ — ^\n{ ^ ^^'' ) as a coefficient in front of the operator G^^,, instead of 
l/Q^. The relevant sum rule takes the form 

1^^(^) + - = ^^Q^T^ + - (') 

(where we once again neglected higher power corrections.) It works well for the 
Euclidean momenta Q^ > 1 GeV^. 

The amplitudes of decays considered above were proportional to ai in the BSW 
language |]I|]. The amplitudes of decays proportional to 02 contain an absolutely 
unknown formfactorfactor, the matrix element < B\lry'^gGauu\'^ >, which cannot be 
detemined using HQET. The sum rules for the decay B^ -^ D^n^ and other decays 
of the type " B ^ D+light meson" in this group will be similar to the above, (with 
the function F{Q'^) instead of l/Q"^) but will include this new formfactor. For the 
decay B -^ J/ihK we have a new function F(Q^) = 2ml \a^2 ^ in the 

sum rule instead of F. 

3 Decay widths 

Let us briefly discuss numerical aspects of our results. We shall concentrate on the 
values of r and the amplitudes ai and a2 that can be directly compared with the 
experimental data. 

For decays B^ -^ D^n^ , B^ —>■ D^p^ we get r ~ —1.5 and ~ —1 respectively. 
For the decays B^ —>■ D*^7i^ , B^ -^ D*^p~ we get r' = r/3 |^. Taking here and 
below Ci ~ 1.12, C2 ~ —0.26, we obtain for these decays ai ~ 1.16, 1.12, 1.08, 1.06 
respectively. 

Consider now other decays using the same method. The discussion below is given 
for orientation only, keeping in mind that the effects unaccounted for in our analysis 
(see section 2) may be important for these decays. If we neglect these effects, we 
obtain for them once again the formulae similar to the one in eq. (3). Consider first 
the decays from the B -^ DD group. Their amplitudes are also proportional to ai and 
can be obtained using the sum rule sketched in section 2. We get for B -^ DD decays 
r ~ —O.9-^0r, where fa is a leptonic decay constant taken to be ~ 170 MeV. We 



immediately see that r ~ —0.8. For B -^ D*D* decays using HQET we obtain r ~ 
— 0-^ 1^2 f 2 ~ —0.16. For B — > D*D decays we get r ~ —0.5. For the corresponding 
Oi factors in the amphtudes we find ai ~ 1.1, 1.04, 1.07 respectively. Note that 
our results for different channels lie between BSW [0] and naive factorization. The 
accuracy of these results is lower than for the previous group of decays, since we 
expect here the perturbative logs and higher corrections can play a bigger role. 

Consider now the decays proportional to the factor 02. Here we shall be extremely 
speculative, since the corresponding analysis is far from being completed. We shall 
only try to indicate what we expect for these decays at the moment, leaving more 
solid statements for the future investigation. The main difficulty here is that we 
do not know the key formfactor— < B\bgGafj.'y'^u\'n' > P°' (and the corresponding 
formfactor with the p-meson). We can try to roughly estimate these formfactors 
from our knowledge of the D meson decays using the symmetry between b and c. 
(Unfortunately, such estimates are very uncertain, though.) Let us completely ignore 
the recoil dependence in the formfactor < B\bgGafj.'y^u\7r > (unlike B^ -^ D^n^ 
there is no justification for that) and parametrize < B\bgGafj.'y^u\7[ > by a number 
m^fj ~ xm'^jj, where x is an unknown constant. Then for the decays of the type 
B^ -^ D°7r° we obtain r ~ —1.6a;, for B^ — »• D*^7r^ we obtain r ~ —0.8a; (the 
difference between the values of r for decays to D and D* is proportional to f'jy/ f^)* 
and we use fo = 170 MeV, /^=220 MeV). The value of x is not known, but the 
experimental data on D seems to indicate that it is below 0.4. If this is indeed the 
case for B decays, then the value of a2 will be suppressed in comparison with the exact 
1/Nc rule for this group of B decays, and can even be equal to zero for sufficiently 
small X. Such a suppression is favored by the recent experimental data |jll|. Future 



calculations are needed to establish x, and at moment we cannot make any definite 
theoretical statement about this group of decays. 

Finally, we note that the same calculation for B — > J /ipK leads to small r due to 
a big leptonic decay constant Fj/^ ~ 300 MeV, 02 ~ 0.12. However in this case there 
exist new difficult problems, due to a large recoil, an enhanced role of higher power 
corrections and higher twists and big continuum contribution ( presumably absent for 
other modes). Moreover, hard gluons can play a significant role in this decay. (The 
sum rule from section 2 has no stability "window" in this case). Thus, we cannot 
exclude the possibility of the rule of discarding l/N^. in this channel yet, neither can 
we confirm it. 

We stress here that the pattern of amplitudes proportional to 02 presented above 
is speculative and is nothing else than an educated guess. A lot of work, especially on 
the determination of chromomagnetic nondiagonal formfactors remains to be done. 



4 Conclusion 

We tried here to draw a general picture for deviations from the naive factorization 
in B decays which stems from nonperturbative QCD. A few aspects requiring further 
clarification are as follows. Higher power corrections must be calculated, perturbative 
logs must be taken into account and the nondiagonal magnetic formfactors must be 
determined. After all this is done the expected accuracy of our results may be 20-30%. 

We also considered K -^ -kti and K — K mixing parameter (see ref. 0), as well 
as inclusive widths of B and D |§, 0, |§] • 

Summarizing, the nonperturbative QCD effects play a key role here in weak 
hadronic decays and can be responsible for the dynamical rule of discarding 1/A^c- 
We now have a general method that allows one to carry out these calculations with 
reasonable accuracy for all decays of B,D and K. 
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